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^ ! Abstract 

The effects of quantum fluctuations due to directional anisotropy and frustration between nearest 

neighbors and next-nearest neighbors of the quantum spin-1/2 Heisenberg antiferromagnet on a 

I . square lattice are investigated using spin-wave expansion. We have calculated the spin-wave energy 

U ■ 
c/3 , dispersion in the entire Brillouin zone, renormalized spin-wave velocities, and the magnetization 

^ ' up to second order in 1/S expansion for the antiferromagnetic Neel and collinear antiferromagnetic 

B' 

stripe phases. It is shown that the second-order corrections become significant with increase in 



frustration. With these corrections magnetizations and spin-wave velocities for both the phases 
become zero at the quantum critical points as expected from other numerical and analytical meth- 



'nI" ' ods. We have shown that the transition between the two ordered phases are always separated by 

> ; 

^ . the disordered paramagnetic phase. 

o 



PACS numbers: 75.10.Jni, 75.40.Mg, 75.50.Ee, 73.43.Nq 



I. INTRODUCTION 

The physics of two-dimensional frustrated spin-1/2 Heisenberg antiferromagnet (HAFM) 
continues to attract considerable attention due to the discovery and availability of new 
magnetic materials such as the layered oxide high-temperature superconductors. ^"^^ These 
systems can be well described by the Heisenberg spin model with nearest neighbor (NN) 
antiferromagnetic coupling Ji and next-nearest neighbor (NNN) antiferromagnetic coupling 
J2. Experimentally by applying high pressures the ground state phase diagram of these 
frustrated spin systems can be explored from low t] = J2/J1 to high rj. For example, 
Li2VOSi04 is an insulating vanadium oxide, with spin s = 1/2 V^"*' ions arranged in square 
lattice planes at the centers of VO4 pyramids. These are linked by Si04 tetrahedra, with Li 
ions occupying the space between the V-0 planes. X-ray diffraction measurements on this 
compound show that the value of 77 decreases by about 40% with increase in pressure from 
zero to 7.6 GPa.^^ Moreover, nuclear magnetic resonance, magnetization, specific heat, and 
muon spin rotation measurements on these compounds (Li2VOSi04, Li2VOGe04, VOM0O4, 
BaCdVO(P04)2) show significant coupling between NN and NNN neighbors.^"® In addition 
these experiments on Li2VOSi04 have shown that it undergoes a phase transition at a low 
temperature (2.8 K) to collinear antiferromagnetic order with magnetic moments lying in 
the a - 6 plane with J2 + ^1 ~ 8.2(1) K and J2/J1 ~ 1.1(1).^'^ 

Quantum spin-1/2 antiferromagnetic Ji — J2 model on a square lattice has been studied 
extensively by various analytical and numerical techniques such as the diagrammatic per- 
turbation theory based on spin-wave expansion^^i^— , modified spin-wave theory^, field the- 
ory2^"— , series expansion^^i"— , exact diagonalization^^, DMRG^"— , effective field theory^i^, 
coupled cluster method^, band-structure calculations,— and Quantum Monte Carlo^"— . It 
is now well known that at low temperatures these systems exhibit new types of magnetic 
order and novel quantum phases.—*^ For J2 = the ground state is antiferromagnetically 
ordered at zero temperature. Addition of next nearest neighbor interactions induces a strong 
frustration and break the antiferromagnetic (AF) order at J2 ~ Ji/2. The competition be- 
tween NN and NNN interactions for the square lattice is characterized by the frustration 
parameter rj = J2/ Ji- It has been found that a disordered paramagnetic phase exists be- 
tween Tjic ~ 0.38 and ri2c ~ 0.60.— 1^ For rj < r/ic the square lattice is AF-ordered whereas 
for rj > ri2c a degenerate collinear antiferromagnetic stripe phase (CAF) emerges. In the 



collinear state the NN spins have a parallel orientation in the vertical direction and antipar- 
allel orientation in the horizontal direction or vice versa. The exact nature of the phase 
transitions and the nature of the intermediate phase are still debatable. It is believed that 
the phase transition from the AF-ordered state to the intermediate paramagnetic state at 
rjic is of second order and from the paramagnetic state to the collinear state at t]2c is of first 
order.— 1^ 

A generalization of the frustrated Ji — J2 model is the Ji — J[ — J2 model where ( = J[/ Ji 
is the directional anisotropy parameter.— i^S It is known that the spatial anisotropy reduces 
the width of the disordered phase. Extensive band structure calculations^ for the vanadium 
phosphate compounds Pb2VO(P04)2, SrZnVO(P04)2, BaZnVO(P04)2, BaCdVO(P04)2 
have shown four different exchange couplings: Ji and J'^^ between the NN and J2 and J2 
between NNN. For example C, ~ 0.7 and J2/ J2 ~ 0.4 were obtained for SrZnVO(P04)2. 
A possible realization of the Ji — J[ — J2 model may be the compound (NO)Cu(N03)3^ 
though recent band-structure calculations show a uniform spin chain model with different 
types of anisotropy and weak interchain couplings^. Within the spin-wave expansion the 
effect of directional anisotropy on the spin-wave energy dispersion and the transverse dy- 
namical structure factor has been studied before.— However, the effect of NNN frustration 
has not been incorporated in that study. 

For the Ji — J[ — J2 model using a higher-order coupled cluster method Bishop et al.— 
reported existence of a quantum triple point (QTP) at C ~ 0.60, 77 ~ 0.33 . Below this 
point they predicted a second-order phase transition between the quantum Neel and stripe 
phases, whereas above it these two phases are separated by an intermediate phase. Existence 
of a QTP has also been reported by other authors^^'^^ where they used effective field theory 
and effective renormalization group approach to obtain a QTP at (^ = 0.51, rj ^ 0.28. In 
a DMRG study it was predicted that there is no intermediate phase (no spin gap) for r) 
lower than 0.287 when C = 1 (isotropic case) .-22. But more recent DMRG calculations have 
concluded that a disordered paramagnetic region persists for all 1] > 0.— 

It should be mentioned that the present Ji — J[ — J2 model was introduced^ as a two- 
dimensional (2D) generalization of the frustrated two-leg ladder. However, the phases of the 
frustrated Ji — J[ — J2 Heisenberg model differs from the phase diagram of the frustrated 
spin-1/2 ladder with rung coupling J[ and diagonal coupling J2.— "^ In case of the frustrated 
spin-1/2 ladder for J[ < 2J2 the ground state is of Haldane type, with two spin-1/2 on the 



rung forming effective spin-1. On the other hand for J[ > 2J2 rung pairs form singlets, 
resulting in the rung-singlet phase.— 

Frustrated two leg ladders share some common features with the present 2D Ji — J( — J2 
model. Using bosonization calculations it has been shown that a spin gap and dimerization 
are also present in this case.^^ However, the presence of this intermediate phase has been 
questioned and a direct transition from the rung singlet to the Haldane phase has been 
reported.— 1^ Yet evidence of such a dimerized intermediate phase in the two leg model was 
found numerically in Refs. 55lJ56| up to a certain value of the interchain coupling. 

One of the main motivations of this work is to investigate (within second-order spin- wave 
expansion) if a disordered paramagnetic region exists for this frustrated, spatially anisotropic 
Ji — J[ — J2 model on a square lattice. We find that the intermediate disordered phase exist 
even for small spatial anisotropics. 

In this work we present a comprehensive study of the effect of zero temperature quantum 
fluctuations on the spin wave energy, spin-wave velocities, and magnetization for the two 
ordered phases of the Ji — J[ — J2 Heisenberg AF on a square lattice. We use spin-wave 
expansion based on Holstein-Primakoff transformation up to second order to numerically 
calculate the physical quantities. Whenever possible we compare our results with available 
experimental data on the systems mentioned above and with other existing analytical or 
numerical results. The paper is organized as follows. Section |TT] provides an introduction to 
the Hamiltonian for the Heisenberg spin-1/2 AF on a spatially anisotropic square lattice. The 
classical ground state configurations of the model and the different phases are then briefly 
discussed. In the next two sections Sec. IIIAI and Sec. IIIBI the spin Hamiltonian is mapped 
to the Hamiltonian of interacting spin-wave excitations (magnons) and spin-wave expansion 
up to second order for spin wave energy, spin-wave velocities, and staggered magnetizations 
are presented for the two ordered phases. These physical quantities for the two phases are 
numerically calculated and the results are plotted and discussed in Section IIIII Finally 
we summarize our results in Section HVl Appendices lAl El and O contain details of the 
formalism. 



II. THE MODEL 



We consider a frustrated S=l/2 antiferromagnet with spatial anisotropy on a Nl x Nl 
square lattice with three types of exchange interactions between spins: Ji along the x (row) 
directions, J[ along the y (column) directions, and J2 along the diagonals. We assume all 
interactions to be antiferromagnetic and positive i.e. Ji, J[, J2 > 0. This Ji — J[ — J2 spin 
system is described by the Heisenberg Hamiltonian 

1 ^ 1 ^ 1 ^ 

«=1 i=l i=l 

where i runs over all lattice sites and i + 6x {Sx = ±1) and i + 6y {6y = ±1) are the nearest 
neighbors to the i-th site along the row and the column direction. The third term represents 
the interaction between the next-nearest neighbors, which are along the diagonals. 

At zero temperature this model exhibits three types of classical ground state configura- 
tions: the Neel state or the (vr, vr) state and the two stripe states which are the columnar 
stripe (vr, 0) and the row stripe (0, vr). The spin orientations of these three states are shown 
in Fig. [H The Neel state breaks the SU(2) and the lattice translational symmetry, but pre- 
serves the fourfold rotational symmetry C4 of the square. The stripe states break SU(2) and 
partial lattice translational symmetries (along one direction). In addition this state breaks 
the invariance under tt/2 real-space rotations C4 to C2. 
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FIG. 1: Classical ground states: (a) AF (7r,7r), (b) Columnar (vTjO), and (c) Row (0,7r) 



The classical ground state energies of these states are determined by treating the spins 



as classical vectors and then minimizing the energy. These are 

1 



Et^/N = --JiS-'z[l + C-H 

EtoiumnJN = -\jiS'z [1 - C + 2r/] , (2) 

EtJN = -\j,S'z[~l + C + H- 

Here C, = J'l/Ji is the directional anisotropy parameter and rj = J2/J1 is the magnetic 
frustration between the NN (row direction) and NNN spins. 2; = 2 is the number of nearest 
neighbor sites. Eq. [2] shows that the classical ground state is either the antiferromagnetic 
Neel (AF) state for t] < (/2 or the columnar antiferromagnetic stripe state (CAF) ii rj > (/2. 
The classical first-order phase transition between the AF and CAF state occurs at the critical 
value 7]^'"'''' = C/2M 

At low temperature quantum fluctuations play a significant role on the phase diagram of 
the system. In the next sections we will consider the classical spins as quantum spins and 
study the role of quantum fluctuations on the AF and CAF ordered phases. We follow a 
standard procedure by first expressing the fluctuations around the "classical" ground state in 
terms of the boson operators using the Holstein-Primakoff transformation.— The quadratic 
term in boson operators corresponds to the linear spin-wave theory, whereas the higher-order 
terms represent spin- wave (magnon) interactions. We keep terms up to second order in 1/S. 
In the next step we calculate the renormalized magnon Green's functions and self-energies. 
Finally, we calculate the magnon energy dispersion, renormalized spin-wave velocities, and 
the staggered magnetization per spin to the leading order in 1/S^ for the AF and CAF 
phases. 

A. AF Phase - Formalism 

For the AF ordered phase NN interactions are between A and B sublattices and NNN 
interactions are between A- A and B-B sublattices. The Hamiltonian in Eq. [1] takes the form: 

i i i 

This Hamiltonian can be mapped into an equivalent Hamiltonian of interacting bosons by 
transforming the spin operators to bosonic creation and annihilation operators a^, a for "up" 



and 6^ , b for "down" sublattices using the Holstein-Priniakoff transformations keeping only 
terms up to the order of l/S"^ 
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Substituting Eqs. |4]we expand the Hamiltonian in powers of 1/S as 



H = -]^N.hS^z{l + C) 



1 



2r7 

T+C 



+ H0 + H1 + H2 + .... 



(5) 



The first term corresponds to the classical energy of the AF ground state (Eq. [2]). Next 
using the spatial Fourier transforms 



N ^ 



-ik-Ri 



e 'flk) Oj = 

k ' " k 



N ^' 



-•^I^Jftk, 



the real space Hamiltonian is transformed to the k-space Hamiltonian. Momentum k is 
defined in the first Brillouin zone (BZ): — tt < k^ < it, — tt < ky < ti (with unit lattice 
spacing). The reduced Brillouin zone contains N/2 k vectors as the unit cell is a magnetic 
supercell consisting of an A-site and a i?-site. 

Furthermore, we diagonalize the quadratic part Ho by transforming the operators Ok and 
6k to magnon operators a^ and /3k using the Bogoliubov (BG) transformations 



«k = ^kttk + "^k/3-k, 6-k = m^al + /k/3-k, 



(6) 



where the coefficients /k and mk are defined as 



- 1 + ek 
• 2ek J 
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■ 2ek J 
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(7) 



with 



7k = 7ik/'^k, 

7ik = [cos(/i;^) + Ccos(A;j^)]/(l + C), 

72k = cos{k:^) cos{ky) , 
2ri 



(8) 



Kk 



1 + C 



(1 -72k)- 



7k is negative in certain parts of the first BZ - so it is essential to keep track of the sign 
of 7k through the function sgn(7k). After these transformations, the quadratic part of the 
Hamiltonian becomes 

Ho = JiSz{l + C) I] /^k (ek - 1) + JiSz{l + 0Y. "^k^k (a]c«k + 4/3k) • (9) 

k k 

The first term is the zero-point energy and the second term represents the excitation energy 
of the magnons within hnear spin- wave theory (LSWT). 

The part Hi corresponds to 1/ S correction to the Hamiltonian. We follow the same pro- 
cedure as described above. The resulting expression after transforming the bosonic operators 
to the magnon operators is 

JiSz{l + ^ J^^ j^^t ^^ ^ ^t^^ j + B^ (^^t ^_^ + (^_^^^ " 

k 

5^ 5g(1 + 2 - 3 - 4)/i/2/3/4 \a\ala^a^vS^^ + Pl^Pl^P^iP^^vi^^ 
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1234 



+ 4al/3U/3-2asVSL + \2a\/3.2asa,vSl + 2/3l,/3_,/3.2asVSL + a\all3i,/3l,vSl 



+ h.c. 



(10) 



In the above equation momenta ki,k2,k3,k4 are abbreviated as 1, 2, 3, and 4. The first 
term in Eq. [10], which is known as the Oguchi correction^ in the literature is obtained by 
setting the products of four boson operators into normal ordered forms with respect to the 
magnon operators, where Ak and Sk are 
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with 
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(13) 
(14) 

(15) 



The second term in Eq. [TU] represents scattering between spin- waves where the delta function 
(5g(1 + 2 — 3 — 4) ensures that momentum is conserved within a reciprocal lattice vector G. 
Explicit forms of the vertex factors V^^]^'^ are given in Appendix lAl 

The second order term, H2 is composed of six boson operators. Before the BG transfor- 
mation H2 is of the following form: 

JiSz{l + 



H, 



{2SyN^ 
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+ i j < 72(2 + 3 — 6)010405020306 + 72(3 — 5 — 6)0^0204030506 

— - (72(3)010402050603 + 72(l)o{o20403050g + /i.e. j +OO6 



(16) 



After transformation to magnon operators a^, /3k the Hamiltonian in normal ordered form 
reduces to 



H2 = '^'^^2sy^^ E [^ik («L«k + /3,^/3k) + C2^ (4/3:^k + /3-k«k) + ... 



(17) 



The dotted terms contribute to higher than second order corrections and are thus omitted 
in our calculations. The coefficients Cik and C2k are given in Appendix [Bl We will find that 
these corrections play a significant role in the magnon energy dispersion and in the phase 
diagram for large frustration and/or small anisotropy. 

The quasiparticle energy E^^ for magnon excitations, measured in units of JiSz{l + C) 
up to second order in 1/S is given as 



E, 



AF 
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(18) 



Expressions for the magnon Green's functions and self-energies are given in Appendix | 



We now define the renormalized spin-wave velocities along the x and y directions at 
the zone boundary using Eq. [18] as \4 = limfc^_^o2Ji5(l + QE^^/k^ with ky = and 
Vy = limfcj^_^o 2 Ji5(l + QE^^/ky with k^ = 0. The renormalization factors are expressed as, 

(19) 

(20) 
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v,y = (C-2r/)-i/^[(C-r;)A + -(1 + 0^2 



(21) 
(22) 

(23) 
(24) 
(25) 

(26) 



The magnetization M defined as the average of the spin operator 5^ on a given sublattice 
(say A) is expressed as 



M Q / t \ c A C _L ^1 ^ ^2 



(25) (25)2' 



where 
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(-E'k + -E'p + -Eq + -Ek+p-q) 

K6) 



k,p,q,[k+p-q] k,p,q,[k+p-q] 



(30) 



E2-(i^p + Eq + Ek+p-q)2 

The zeroth-order term AS corresponds to the reduction of magnetization within LSWT, Mi 
term corresponds to the first-order 1/S correction, and M2 is the second-order correction. 
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B. CAF Phase - Formalism 



1. Hamiltonian 



In the CAF phase "up" and "down" spins interact along the row directions (NN couphng) 
and also along the diagonals (NNN coupling) whereas "up" -"up" and "down" -"down" spins 
interact along the column direction (NN coupling). The Hamiltonian for this phase is de- 
scribed by 



H - J^Z^^i ■ Sj+5, + 2 "^^ ^ [^^ ' ^^+^'^+^y "^ ^' ' ^i+Sx+Sy 



J2Y.St-Sls,. (31) 



The Hamiltonians for the AF and the CAF ordered phases (Eq. [3] and Eq. [3T]) show the 
similarity between these two phases. In the AF-phase J2 interactions play the role of J[ 
interactions in the CAF phase. For the CAF-phase the structure factors 7jj^, 725^. along with 
other quantities required for the calculations are defined as 
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The coefficients for the Oguchi correction that appear in the Hamiltonian Hi are 
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Hq, Hi, and H2 can be expressed in the same forms as in Eqs. [9l [TO], and [T7| with the 
new coefficients A'^^, B^, C[y^, Cg^ and with the replacement ( ^ 2rj. The expressions for the 
two vertex factors V'^^\ V'^^^ and the coefficients C(j^, Cgi^ are given in Appendix O As an 
example for the CAF phase Eq. IH] takes the form: 



Ho = JiSz{l + 2r7) Y, 4 (4 - 1) + JiSz{l + 2r7) J^ n'^e'^ (aj^ak + PIP^ 



(38) 



The quasiparticle energy E^^^ for magnon excitations, measured in units of 715*^(1 + 2?]) 
up to second order in 1/ S is given as 



E, 
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(39) 



The renormalized spin-wave velocities along the x and y directions for this phase are 
defined as V^ = limfc^^o2JiS'(l + 2ri)E^^^ /k^ with ky = and Vy = \im.ky^o2JiS{l + 
2r))E^^^/ky with k^ = 0. The renormalization factors are, 
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2Ei 



(40) 
(41) 

(42) 
(43) 
(44) 
(45) 
(46) 

(47) 



III. RESULTS 



AF Phase 



1. Spin-wave energy dispersion 

We numerically evaluate Eq. [18] to obtain the spin-wave energy 2715(1 + QE^^ as a 
function of momentum for several values of C, and tj. For the numerical summation, the first 
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BZ is divided into N^ meshes with A^^ = 64 and then 4096 points of p and 4096 points of 
q are summed up to evaluate the third term in Eq. [181 For some of the cases we have used 
Nl = 96 for better accuracies. 

Figure 12] shows a comparison between the results from LSWT (long-dashed lines), first- 
order (dotted lines) and second-order corrections (solid lines) to the spin-wave energy. 
Fig. |2]^a) shows the spin-wave energies for isotropic coupling {( = 1) with 77 = and 77 = 0.3 
and Fig. Et^b) shows the results with two different values of anisotropy parameter ( = 1 and 
0.4 for ?7 = 0. We find that in the entire BZ both the first (dotted lines) and second order 
corrections (solid lines) make the spin-wave energy larger and the corrections from LSWT 
(long-dashed lines) are significant for all cases. 

In Fig. iniwe show the spin- wave energy results with second-order corrections for different 
values of ( and rj. The spin- wave energy curve for the isotropic coupling (^ = 1 with 
1] = was reported earlier.— The dispersion along (vr/2, 7r/2)-(7r, 0) is flat within LSWT and 
1/S correction (See Fig. |2]for example). The second-order corrections make the excitation 
energies at (7r,0) smaller than the energies at (7r/2,7r/2). Our results for spin- wave energy 
with frustration and with anisotropic couplings are new. The dip in the magnon energy at 
(tt, 0) increases with increase in frustration rj. Experimentally this can provide a measure of 
the strength of NNN frustration. 

In Fig. m we show the effect of the directional anisotropy parameter ( on the spin- wave 
energy (with second-order corrections). Similar to Fig. |3]we find that the dip in the energy 
at (vr, 0) increasing values of (. 

Recently using neutron scattering measurements on copper deuteroformate tetradeurate 
(CFTD), a real two dimensional Heisenberg AF with weak interplane interactions (~ 10^^ — 
10~^Ji) magnon energies have been obtained for the entire BZ.-i^. It was found that the 
energies at (vr, 0) is 13.5180 meV (with estimated error of 0.1641 meV), which is about 7(1)% 
smaller than the energy 14.4880 meV (with estimated error of 0.0647 meV) at (vr/2, vr/2).-'^ 
The coupling Ji is estimated to be 6.19 meV. This local minimum at (vr, 0) is due to quantum 
fluctuations and may be due to multimagnon processes (entanglement of spins on neighboring 
sites) at this zone boundary. Series expansion around the Ising limit'^^ and Quantum Monte 
Carlo methods^ have accounted for all of the experimental data. But these numerical 
methods do not provide any insight into the physics at this zone boundary. To test our 
numerical procedure we systematically calculate the values of E^^^JJi and E^J/2-k/2)I --^"^ ^^^ 
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FIG. 2: Spin-wave energy Eq^ / Ji obtained from LSWT (long-dashed lines), with first-order (dot- 
ted lines) and second-order corrections (solid lines) for the AF-ordered phase. Figure (a) is for 
isotropic coupling C, = 1 with r/ = and r] = 0.3 and Fig. (b) is for r/ = (no frustration) and 
with two different values of anisotropy parameter C = 1 ^'^^ 0-4. Both the first and second order 
corrections make the spin-wave energy larger and the corrections from LSWT (long-dashed lines) 
are significant in the entire BZ. (color online) 

Nl = 36, 48, 64, 96 and 128. The convergence of our results are very good as shown in FigEi 
We extrapolate these results using the fitting function A + B/Nl + C/N^ + D/N^ to obtain 
A for Nl — i- c>o and reproduce the numerical results Ef^^JJi ^ 2.3585, E^J/2-k/2)I ■-^'^ ~ 
2.3908 reported earlier— with a 1.4% decrease between these two energy values. Recently 
a third-order in 1/S expansion has been done to obtain the spectrum of short-wavelength 
magnons^^*^ where it was shown that the 1/S series converges slowly near the wave- vector 
(7r,0). With the third order correction the excitation energy at (vr,0) was found to be 
3.2% smaller than at (7r/2,7r/2). This result for the energy difference still falls short of 
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FIG. 3: Spinwave energy EQ^ / Ji for the AF ordered phase with second-order correctfons is plotted 
for different values of C, and rj. The dispersion along (vr/2, 7r/2)-(7r, 0) is flat within LSWT and 
1/S correction. The second-order corrections make the excitation energies at (vr, 0) smaller than 
the energies at ('7r/2, vr/2) for all cases. With increase in NNN frustration r] (for a fixed value of 
the directional parameter C,) the dip in the magnon energy at (vr, 0) increases. This can provide a 
measure of the strength of NNN frustration, (color online) 

the experimental result of 7%. This suggests that the inclusion of correction to even third 
order in 1/S is insufficient to explain this energy difference. It should be noted that other 
interactions e.g. ring exchange interactions have been proposed to play a role in these 
compounda^i^i^. 

It may be interesting to study the effects of small NNN frustration and small anisotropy 
on the energies at these two zone boundaries. Table [T] shows our extrapolated values of 
E^^qJ Ji and -S(^/2 7r/2)/'^i ^^'^ ^^^ percentage changes for small frustrations rj = 0.01,0.02 
and a small directional anisotropy ( = 0.98. Calculations are done with lattice sizes N^ = 
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FIG. 4: Effect of directional anisotropy parameter (^ on the spin- wave energy. The calculated 
spin- wave energy is with second-order corrections. Similar to Fig. [3] the dip in the energy at (tt, 0) 
increases with increase in the values of (. (color online) 
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FIG. 5: Convergence of E^J^^-^/Ji and E^J/2,^/2)/Ji is shown with 1/Nl for Nl = 48,64,96 and 
128. We extrapolate these results using the fitting function A + B/Nl + C /Nf^ + D/N^ to obtain 
^KO)/-^! = 2-3585, ^(t^^^/a)/-/! = 2.3908 for Nl ^ oo. 

48, 64, 96 and 128 and the results are extrapolated to N^ — )■ oo using the fitting function 
A + B/Nl + C/Nl + D/Nl. We show that a small frustration (for example rj = 0.02) for 
the isotropic coupling causes a noticeable difference (2.8% within second-order spin-wave 
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expansion) in energies between E^J^J Ji and E^J/2-k/2)I ■-^'^- These features can be explored 
experimentally using neutron scattering measurements with compounds that can be modeled 
by the Ji — J( — J2 Heisenberg antiferromagnet. 

TABLE I: Energies E^^^JJi and Ef^J/2-ni2)l '^'^ ^^'^ different values of C, and r/. 



A''l -^ 00 




EtI,)lJ^ 


^(t/2,V2)/'^l 


^E'^'/Etl,,,.,.) 




LSWT 


2.0000 


2.0000 


0% 


C=l,7? = 


1/S 


2.3159 


2.3159 


0% 




1/S2 


2.3585 


2.3908 


1.4% 




LSWT 


1.9600 


1.9800 


1.0% 


^=1,7^ = 0.01 


1/S 


2.2800 


2.2980 


0.8% 




1/S2 


2.3221 


2.3753 


2.2% 




LSWT 


1.9200 


1.9600 


2.0% 


( = 1,7^ = 0.02 


1/S 


2.2443 


2.2801 


1.6% 




1/S2 


2.2886 


2.3536 


2.8% 




LSWT 


1.9799 


1.9800 


0% 


( = 0.98,r? = 


1/S 


2.2926 


2.2928 


0% 




1/S2 


2.3348 


2.3680 


1.4% 



2. Renormalized spin-wave velocities 



We calculate the spin- wave velocity renormalization factors Z^J" , Z^ along the x and y 
directions from Eq. [201 For the second-order correction terms V2xi V2y, we consider lattice size 
Nl = 72 and evaluate [Saci(k, Ek) — je"]/^^ with A;^ = 71 /N^. V2y is obtained similarly. For 
the isotropic case C = 1 and with 77 = we find the second-order correction V2x = V2y = 0.021 
which is in excellent agreement with results reported earlier.— "i^i^i The results from our 
calculations with increase in rj are shown in Fig. El We find that the velocities steadily 
decrease with increase in frustration and finally becomes zero close to the quantum critical 
points ?7ic for the AF-phase. Second-order corrections are significant to stabilize the velocities 
as with first-order corrections these velocities diverge with increase in frustration (similar 
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to the case with magnetization discussed later). We also notice that the difference between 
the renormalization factors Z!^^ and Z^^ diminishes with increase in frustration. 
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FIG. 6: Renormalization factors Z:^^ and Z:^^ for the spin- wave velocities are plotted with frus- 
tration r] for different values of C. The velocities steadily decrease with increase in NNN frustration 
and finally becomes zero close to the quantum critical transition points rjic for the AF-phase. (color 
online) 



3. Staggered Magnetization 



We obtain the staggered magnetization Map for the AF phase with several values of ( 
and T] from Eq. [27] by numerically evaluating Eqs. [25H5U1 Especially to obtain the second 
order correction term M2 we sum up the values of Nf^/A points of k in the 1/4-part of the 
first BZ and Nl points of p and q in the first BZ, with Nl = 36 lattice sites (total of about 
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544.2 million points for each ( and 77). Except for small spatial anisotropy (, M2 values 
start from a positive small number and then switch sign and become negative with increase 
in frustration rj. However, for small (, say ( = 0.2 M2 starts from a small negative number 
(~ —0.005) and remains negative with increase in rj. Figure [7| shows the magnetization with 
increase in the frustration parameter rj = J2/J1 for several values of the spatial anisotropy 
parameter ( = J[/Ji = 0.2,0.4,0.6,0.8, and 1.0. For each ( three different curves are 
plotted: the long-dashed lines represent LSWT prediction, the dotted lines include the 
first-order (l/S) correction to the LSWT results, and the solid lines represent corrections 
up to second-order [l/S"^) to the LSWT results. With increase in frustration the dotted 
curves diverge. However, 1/S^ corrections (M2) significantly increase with frustration and 
stabilize the apparent divergence of the magnetization. We find that the magnetization with 
second-order corrections decreases steadily at first and then sharply drops to zero. As an 
example, for the isotropic case (C = ^,V = 0); Maf starts from 0.307 and then decreases 
till rj ^ 0.32 and finally becomes zero at the critical point rjic ~ 0.41. For this case we 



reproduce the magnetization plot obtained in Ref. Il3|. Other values of ( show the same 
trend except for small (. For ( = 0.2 we find Map steadily decreases from 0.21 at rj = 
to 0.19 at r/ ^ 0.054 and then slightly increases to 0.195 at r^ ^ 0.068. Finally it sharply 
drops to zero at t]ic ~ 0.084. This feature has not been observed before and may be an 
artifact of the spin-wave expansion showing the limitation of this method for small ( (the 
system becomes essentially one- dimensional as C — ^ 0). It may be interesting to verify this 
by series expansion or other analytical or numerical methods. Note that for all cases second- 
order corrections increase the critical value of rjic from the LSWT predictions. Our values of 
magnetization for 77 = (no NNN frustration), M^p = 0.307 agrees with previously obtained 
values from spin-wave expansioni^"— 1^, series expansion^"—, and experimental results for 
K2NiF4, K2MnF4, Rb2MnF4 and other system8^'S~-ii. 

The ground-state magnetization per spin is reduced from its classical value S = 1/2 by 
zero-point quantum fluctuations. This "spin reduction" Aaf = 0.5 — M^p is plotted for 
different values of C for 77 = in Fig. [HI Second-order corrections (solid line) change the 
values of M^p slightly from the LSWT predictions (dashed line). The fluctuations increase 
with decreasing values of C, suggesting that spin-wave expansion for S = 1/2 is not applicable 
for (^ < 0.1 as the system essentially becomes one dimensional. 

In the inset of Fig. [8] we show the spin deviation with l/( for (^ = 0.1 to 1. We find 
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FIG. 7: Staggered magnetization M\p is shown for the AF ordered phase with frustration i] = J2/J1 
with different values of spatial anisotropy ( = J{/Ji. For each (^ results from LSWT (long-dashed 
lines), with first-order (dotted lines) and with second-order corrections (solid lines) are plotted. 
With increase in rj the dotted curves diverge. Second-order 1/S^ corrections become significant 
for large rj and they stabilize the apparent divergence of the magnetization. Magnetizations with 
1/S^ corrections decrease steadily and then sharply drop to zero. For example in the isotropic case 
i.e. C = 1) -^AF starts from 0.307 and then decreases till rj ~ 0.32 and finally becomes zero at the 
critical point r/ic ~ 0.41. However for small Ci say (" = 0.2 we find M\p to steadily decrease from 
0.21 at ?7 = to 0.19 at r] ^ 0.054 and then slightly increases to 0.195 at rj ^ 0.068. Finally it 
sharply drops to zero at rjic ~ 0.084. Note that for all cases second-order corrections increase the 
critical value of r/ic from the LSWT predictions, (color online) 

that the function /(C) = 0.16 + 0.029C"^ - 0.00079C^^ is a good representation of Aaf for 
this range of (. The fitted curve (dashed hne) is shown in the inset along with the actual 
numerical results (solid line) with the second-order corrections. 
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FIG. 8: Spin deviation Aaf = 0.5— M^p from the classical value of 0.5 is plotted for the AF ordered 
phase (with no NNN interaction, i.e. rj = 0) for different values of spatial anisotropy C. Dashed 
line is LSWT prediction whereas the solid line includes 1/S'^ corrections to LSWT results. The 
fluctuation increases with decreasing values of C, suggesting that spin-wave expansion for S = 1/2 
is unreliable for ^ < 0.1. In the inset we show the fluctuations with l/(^ for ( = 0.1 to 1. The 
function /(C) = 0.16 + 0.029C~^ - 0.00079C~^ is a good representation of Aaf for ( = 0.1 - 1.0. 
The fitted curve is shown in the inset along with the actual numerical results with the second-order 
corrections, (color online) 

B. CAF Phase 

1. Spin-wave energy dispersion 

We numerically evaluate Eq. [39] with A^^ = 72 lattice size to obtain the spin- wave energy 
2 JiS'(l + 2?7)i^^^^ as a function of momentum for several values of ( and rj. The calculations 
are similar to the AF-phase. Figure M shows the spin-wave energy E^^^ / Ji with second- 
order corrections (solid lines) for the entire Brillouin zone of the CAF-ordered phase. Two 
different values of directional anisotropy C, = 0.9 and 0.4 for NNN frustration t] = 1 are 
chosen. Results obtained from linear spin-wave theory (long-dashed lines) and with only 
first-order (dotted lines) are also shown for comparison. 1/S and 1/S^ corrections increase 
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FIG. 9: Spin-wave energy E^ /Ji results obtained from linear spin-wave theory (long-dashed 
lines), with first-order (dotted lines), and second-order corrections (solid lines) for the entire Bril- 
louin zone of the CAF-ordered phase. For fixed NNN frustration 77 = 1 two different values of 
C = 0.9 and 0.4 are chosen. 1/S and l/S'^ corrections increase the spin-wave energy of the or- 
dered phase from the linear-spin wave theory results. For (^ = 0.4 second order corrections are 
insignificant compared to the first order 1/S corrections. However, for (" = 0.9 1/S^ corrections 
lower the spin-wave energy from the first-order corrections. Spin-wave energy shows three peaks, 
the maximum being at (7r/2,0). The second small peak is at (0.5147r, 0.4867r) and the third peak 
occurs at (7r,7r/2). (color online) 

the energy of the ordered phase from the LSWT results. We find that the second-order 
corrections to the magnon energy are not significant from the energy obtained with first- 
order corrections for small (. However for large (, say ( = 0.9 1/S^ corrections lower the 
spin- wave energy from the first-order 1/S corrections. 

Figure [10] shows the effect of frustration 77 for a fixed value of spatial anisotropy ( = 
0.6. Second-order corrections are negligible compared to the first-order corrections, which 
significantly enhance the LSWT results. In both Fig. [9] and Fig. [10] the spin-wave energy 
vanishes at the wave- vector (vr, 0) as expected for the CAF phase. We find three peaks in 
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the magnon energy, the maximum being at (vr/2, 0). The second small peak in energy is at 
(0.5147r, 0.4867r) and the third peak occurs at (tt, 7r/2). 
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FIG. 10: Effect of frustration r/ on the spin-wave energy in the CAF-phase. Second-order correc- 
tions are negligible compared to the first-order corrections. However, 1/S corrections significantly 
enhance the spin-wave energy obtained from LSWT results, (color online) 



2. Renormalized velocities 



Renormalization factors Z^J^^, Z^^^ along the x and y directions are obtained from Eq. 
with second-order corrections. The results are shown in Fig. [TT] As we expect similar to 
the AF-phase the velocities steadily decrease with increase in rj and finally becomes zero 
close to the quantum transition points 772c for the CAF-phase. Second-order corrections are 
significant to stabilize the velocities as with first-order corrections these velocities diverge 
with increase in frustration (similar to the case with magnetization discussed later). For 
small C (C = 0.2, 0.4 in figure) we find that Z^^^ slightly increases and then sharply drops 
to zero. We increase the lattice size to N^ = 96 to check the accuracy of our calculation. 
We find no changes in our plot. It may be interesting to verify this with series expansion 
or other analytical or numerical methods. Our numerical method based on the spin-wave 
expansion for the CAF phase is not reliable for ( > 0.95 - so we have not been able to obtain 
the renormalized spin-wave velocities for the case with C = 1 and rj = 1 (more discussed in 
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FIG. 11: Renormalization factors -^^ and Z'^ are plotted with frustration t] for different values 
of C- The velocities steadily decrease with increase in NNN frustration and finally become vanish 
close to the quantum critical transition points 7/2c for the CAF-phase. Numerical calculations are 
done with lattice size Ni^ = 72. For C, = 0.2,0.4 we find that Z^ slightly increases and then 
sharply drops to zero. Increasing the lattice size to A''^ = 96 does not change our results, (color 
online) 



3. Staggered Magnetization 

Similar to the AF-phase the staggered magnetization Mcaf for the CAF phase with sev- 
eral values of ( and r] are obtained by summing over points in the first BZ with Nl = 36 
lattice sites. Except for large spatial anisotropy (, M2 values start from a small positive 
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number and then switches sign and become negative with increase in frustration 77. However, 
for large (, say for C = 0.8 M2 corrections are always negative. Figure [12] shows the magneti- 
zation with increase in frustration parameter r] for several values oi ( = 0.2, 0.4, 0.6, 0.8, 0.95. 
For each ( three different curves are plotted: LSWT results (long-dashed line), first-order 
corrections (dotted line), and second-order corrections (solid line) to the LSWT results. 
Similar to the AF-phase the dotted curves diverge with increase in frustration. However, 
1/S^ corrections (M2) significantly increase with frustration and stabilize the magnetization 
and finally make it zero. We find that Mcaf decreases steadily at first and then sharply 
drops to zero. As an example, for the C = 0.2 Mcaf starts from 0.371 and then decreases 
till rj ^ 0.12 and sharply drops to zero at the critical point r]2c ~ 0.116. With increase in ^ 
the values of the critical points ri2c differ more from LSWT predictions. 

We also find that starting from ( = 0.95 the spin deviation Acaf = 0.5 — M^^p increases 
substantially as we approach the isotropic limit ( = 1. This is shown in Fig. [131 Acaf 
from the LSWT theory remains smooth (dashed lines in Fig. [T3l) . Both the first (Mi) and 
second order (M2) corrections increase rapidly for C > 0.95. This increase is due to the 
fact that e|^ — 7- as C — ^ 1- We have not found a numerical way to regulate it. Instead we 
used extrapolation to obtain values of Acaf beyond ( = 0.95. Inset of Fig. [13] shows both 
the exact data (solid line) and the extrapolated curve (dashed line). With the extrapolated 
curve we obtain Acaf ~ 0.20 for ^ = 1, which gives Mcaf ~ 0.30 for the isotropic limit. 
This is in good agreement with the recent neutron scattering measurements data of the order 
parameter Mcaf = 0.31(2) for Li2VOSi04, which is believed to be a S' = 1/2 frustrated 
antiferromagnet on a square lattice with J2 ~ Ji.-~- 

C. AF and CAF ordered phases - Phase Diagram 

Staggered magnetizations of a spatially anisotropic frustrated spin- 1/2 Heisenberg anti- 
ferromagnet on a square lattice is presented in Fig. [H] for both AF and CAF ordered phases. 
We find the staggered magnetization M^p for r^ = to decrease with decrease in anisotropy 
( in the AF-phase. On the other hand, M^^^p increases with decrease in ( for r] = 1. Our 
results for M^p = 0.307 for the AF-ordered phase and M^^p = 0.30 for the CAF-ordered 
phase are in excellent agreement with existing experimental data on these systems. Further- 
more, we find that in both the phases the second order corrections play a significant role to 
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FIG. 12: Staggered magnetization Mcaf for the CAF ordered phase is plotted with rj for different 
values of C- For each value of C three different curves are shown: long-dashed line is the prediction 
from LSWT, dotted line is the first-order correction, and the solid line includes corrections up to 
second-order. In all cases first-order corrections diverge for some value of r/. However, second- 
order (1/S^) corrections become significant and stabilize the magnetization. Similar to the AF- 
phase Mcaf with second order corrections decreases steadily and then sharply drops to zero. For 
example with C = 0.2, Mcaf starts from 0.371 at r/ = 1 and then decreases till r] ~ 0.12 and sharply 
drops to zero at the critical point rj2c ~ 0.116. For ( more than 0.95 the fluctuations become too 
large (see Fig. llSp - in that case our spin- wave expansion becomes invalid (see text). With increase 
in ( the values of the critical points r]2c differ more from the LSWT predictions. (color online) 

stabilize the magnetization. Staggered magnetizations become zero in both the phases at 
the critical values rjic and ri2c for each value of (. 

Phase diagram for the Ji — J( — J2 model is displayed in Fig. [151 The solid lines indicate 
the critical points rjic for the AF and ri2c for the CAF phases. The dashed line is the 
classical first-order phase transition line between the two phases. Our spin-wave expansion 
for the CAF phase becomes unreliable for 77 > 0.95. Thus we extrapolate our data to obtain 
V2c = 0.58 for ( = 1. The dotted line in Fig. [15] is the extrapolated curve. This result 
r]2c = 0.58 is in good agreement with the expected value of ~ 0.60 for the isotropic case. 
Figure [15] shows that the spin-gap (772c — f]ic) increases with increase in (. 
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FIG. 13: Spin deviation Acaf is plotted for the CAF ordered phase {rj = 1). Dashed hne is LSWT 
results and the solid line is with the l/S*^ corrections. For (" more than 0.95 spin-wave expansion 
becomes unreliable as Acaf as the first and second-order corrections to Acaf increase rapidly. 
Thus with the 1/S^ corrections we extrapolated the values to obtain Acaf ~ 0.20 for C = 1- This 
gives McAF ~ 0.30 for the isotropic limit. This is in good agreement with existing experimental 
results (see text). Inset shows both the exact data (solid line) and the extrapolated curve (dashed 
line), (color online) 

Within our spin-wave expansion we do not find any quantum triple point for any values 
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and 



44 . Instead from our cal- 



of C and rj. This is in contrary to the findings in Refs. 
culations we find that there are two ordered phases separated by the magnetically disordered 
phase. Our proposed phase diagram is consistent with the phase diagram obtained by the 
DMRG calculations^, exact diagonalization method^, and the results from the continuum 
limit of the present modeP^. Our results are also in accord with the numerical evidence of 
a dimerized intermediate phase in the frustrated two leg model up to a certain value of the 
interchain coupling.— i^ 
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FIG. 14: Staggered magnetization M is plotted for both the AF and CAF ordered phase (with 
second-order corrections). For the AF phase the different values of (" are 0.2, 0.4, 0.6, 0.8, and 1.0 
and for the CAF phase the values are 0.2, 0.4, 0.6, 0.8, and 0.95. Our numerical approach using 
spin- wave expansion is not reliable for the CAF phase for C larger than 0.95 (see text for details). 
For both the phases M become zero at some critical values of the NNN frustration parameter ij. 
We also find that the spin-gap increases with increase in r]. (color online) 

IV. CONCLUSIONS 

In this work for an antiferromagnetic square lattice we have provided a comprehensive 
study of the effects of quantum fluctuations due to spatial anisotropy and frustration be- 
tween nearest and next-nearest neighbors on the low-temperature thermodynamic properties 
of the two ordered phases of the system. Using second-order spin-wave expansion we have 
calculated the spin wave energy in the entire Brillouin zone, renormalized spin-wave veloci- 
ties, and the magnetizations for the antiferromagnetic Neel and columnar antiferromagnetic 
phases. We have found that the second-order corrections contribute significantly to stabilize 
the quantum phase diagram of the system as frustration between the spins increase. As 
expected from linear spin wave theory magnetization becomes zero at the quantum critical 
points. However, the second-order corrections slightly extend the region of the AF-order. 

Our results for the spin-wave energies are compared with the recent experimental re- 
sults using neutron scattering for CFTD.- With our second-order spin- wave expansion we 
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FIG. 15: Phase diagram for the Ji — J( — J^ model. The sohd lines indicate the critical points 
r/ic for the AF and r/2c for the CAF phases. Our spin- wave expansion for the CAF phase becomes 
unreliable for 77 > 0.95. We extrapolate our data to find r\2c = 0.58 for ( = I. The dotted line 
is the extrapolated curve. The dashed line in the middle represents the classical first-order phase 
transition line 77^ = C/2. The spin-gap {r]2c — rjic) increases with the anisotropy parameter ^. 
(color online) 

have reproduced the previous numerical resuhs that the spin-wave energy at (vr, 0) is about 
1.4% smaller that at (7r/2,7r/2). This result falls short of the experimental result. Further- 
more, we find that the dip in spin wave energy at (tt, 0) increases with increase in NNN 
frustration. This can provide a measure of the effect of frustration experimentally. For a 
few values of small frustration and anisotropy we have explicitly calculated the percentage 
changes between the spin-wave energies at {tt,0) and at {tt/2,tt/2). We have shown how 
the renormalized spin-wave velocities along the row and column direction change with frus- 
tration. Both these velocities become zero close to the critical transition points. For the 
AF-ordered phase we have also calculated the spin deviation from the classical value of 0.5 
with no NNN coupling for different values of directional anisotropies and have obtained an 
empirical equation based on our numerical data. 

For the CAF-ordered phase we have obtained similar results. The magnetization becomes 
zero at the quantum critical points as frustration increases. For ( < 0.95 our calculations 
produce correct results but our present numerical approach is not reliable for ( > 0.95. 
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Thus we were not able to find the thermodynamic properties ioi ( = 1. Based on our data 
we have extrapolated the magnetization for the case ( = 1,1] = 1 and found it to be 0.30 
which is in good agreement with existing experimental result^"-. Our extrapolated value 
of the quantum critical point ri2c = 0.58 for C = 1,?7 = 1 is also in good agreement with 
the expected value ~ 0.60. We have not found much experimental data on this system to 
compare with our other results such as the spin-wave energy dispersion in the entire BZ and 
the spin-wave velocities. 

Finally we combined our results for the magnetization of the two phases with different 
directional anisotropics to obtain the complete magnetic phase diagram of the system. We 
have found that two ordered phases are always separated by the disordered paramagnetic 
phase. Our proposed phase diagram is consistent with the phase diagram obtained from the 
DMRG calculations^, exact diagonalization method^, and the results from the continuum 
limit of the present model^^. Our results are also in accord with the numerical evidence 
of a dimerized intermediate phase in the frustrated two leg model up to a certain value 
of the interchain coupling. In summary with our present approach based on second-order 
spin-wave expansion we do not find existence of quantum triple points for any values of C 
and rj. 
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Appendix A: Vertex factors for the AF phase 



The six vertex factors for the AF-phase are given below. 



V, 



(1) 



1234 



+ 



7i(l - A)xiXi + 7i(l - 3)a:ia;3 + 7i(2 - A)x2Xi + 7i(2 - 3)a;2a;3 
- 7i(l)a;i + 7i(2)a;2 + 7i(3)a;3 + 71(4)0:4 + 71(2 - 3 - 4)0:2X3X4 

7i(l — 3 — 4)a;ia;3a;4 + 71 (4 — 2 — 1)0:1X2X4 + 71 (3 — 2 — 1)^1X2X3 
2ri 



1 + C 



/1234 1 + Sgn(7G)xiX2X3X4 



V, 



(2) 



1234 



7i(2 - 4)a;ia;3 + 71(1 - 4)^2X3 + 71(2 - 3)xia;4 + 71(1 - 3)x2X4 

- 71(2)^1X3X4 + 71(1)2^2X3X4 + 7l(4)xiX2X3 + 7i(3)xiX2X4 

7i(2 - 3 - 4)a;i + 71(1 - 3 - 4)a;2 + 71 (4 - 2 - l)a;3 + 7i(3 - 2 - l)x4 
2r] 



- (^Y— :)/i234 X1X2X3X4 + sgn(7G) 



V, 



(3) 



1234 



7i(2 - 4) + 7i(l - 3)xiX2X3X4 + 71(1 - 4)a;iX2 + 7i(2 - 3)x3X4 

1 r 

- 7i(2)x4 + 71(1)3:1X2X4 + 7i(2 - 3 - 4)a;3 + 7i(l - 3 - 4)^1X2X3 



7i(4)a;2 + 7i(3)x2X3X4 + 71 (4 - 2 - l)a;i + 71(3 - 2 - 1)0:1X3X4 
2r] 



1 + C 



/: 



1234 



X2X4 + Sgn(7G)xiX3 



V, 



(4) 



1234 



+ 



-71(2 - 4)a;4 - 7i(l - 4)0:1X2X4 - 7i(2 - 3)x3 - 7i(l - 3)xiX2X3 
- 7i(2) + 71(1)^1X2 + 7i(3)x2X3 + 7i(4)x2X4 

7i(2 - 3 - 4)0:3X4 + 7i(l - 3 - 4)0:1X2X3X4 + 7i(3 - 2 - l)xia;3 + 71 (4 - 2 
2r] 



(Al) 



(A2) 



(A3) 



1 + C 



/1234 X2 + Sgn(7G)xiX3X4 



l)XiX4 

(A4) 



V, 



(5) 



1234 



+ 
+ 



-71(2 - 4)a;i - 7i(2 - 3)xiX3X4 - 7i(l - 4)a;2 - 7l018i(l - 3)x2X3X4 

- 71(2)^1X4 + 71(1)0:2X4 + 7i(4)xiX2 + 7i(3)XiX2X3X4 

7i(2 - 3 - 4)a;ia;3 + 7i(l - 3 - 4)^2X3 + 71 (4 - 2 - 1) + 71 (3 - 2 - 1)0:3X4 
2?7 



1 + C 



/: 



1234 



X1X2X4 + Sgn(7G)x3 



(A5) 



V, 



(6) 



1234 



+ 



7i(2 - 4)^2X3 + 7i(2 - 3)x2X4 + 7i(l - 3)xia;4 + 7i(l - 4)a;iX3 
- 71(2)^2X3X4 + 7i(3)x4 + 7i(2 - 3 - 4)x2 + 71 (3 - 2 - 1)^1X2X4 



7i(l)xiX3X4 + 7i(4)x3 + 7i(l - 3 - 4)xi + 71 (4 - 2 - l)xiX2X3 
2r] 



c 



/1234 X3X4 + sgn(7G)xiX2 



(A6) 
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with 



/: 



1234 



72(1 -3) +72(1 -4) +72(2 -3) +72(2 -4) -72(1) -72(2) -72(3) -72(4) . (A7) 



Appendix B: Green's function and Magnon Self energy for the AF phase 

The time-ordered magnon Green's functions are defined as 

G,„(k,t) = -z(TK(i)4(0))), G'^/3(k,t) = -z(T(/3lk(t)/3_k(0))), 

Considering Hq as the unperturbed Hamiltonian the Fourier transformed unperturbed prop- 
agators are given as 



GL(k,^) = [uj-E, + t6]-\ GUk,u) = [-u-E, + z6]'' 



Gl^iKu) = G%ik,u)=0, 



(Bl) 
(B2) 



with 6 — )■ 0-I-. The magnon energy E^ = KkCk is measured in units of JiSz{l + (). The 
graphical representation of the Green functions are shown in Fig.fTBlfa). The full propagators 



-I- 



(a) 



a 



> 



a 



(^ 



-<- 



P 



(b) 



q 



■<- 




k+p-q \ k 

—< ) < + 



(4)" --.>■-' (4) 
V P V 




.^6^^^7<-'V^^ 



V 



(c) 



q 



■^ 




k+p-q \ ,, 
< ) <--■ + 



^(4r^^,>.'^(6) 




v^'^"^\<-'\^'^ 



FIG. 16: (a) The solid and the dashed lines correspond to the a and /? propagators. Second-order 
diagrams for the self-energies SQ,Q,(k,c;j) and SQ,^(k, c^j) are shown in (b) and (c). V^"^' ,V^^',V^^' 
are the vertex factors (see text). 
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Gjj(k, w) satisfy the matrix Dyson equation: 

G., (k, io) = G°,.(k, i0) + Y, GLiK Uj)^mn{K Oo)GnAK 

•mn 

where the self-energy Sjj(k can be expressed in powers of 1/(25') as 



T.ij(k,uj) 



(2), 



{2S) '^ ^^'^^^ (25)2"*^ 



k,u) + 



(B3) 



(B4) 



Using the relations 



.(6) 



^k+p-q],q,p,k ~ ^ShITg) V|^^p^q [k+p_q] , ^q,[k+p-q],k,p ~ SSn(7G)^k,p,q,[k+p-q] ' 

the first and second order self-energies are written as 



Si^i(k,c.) = : 


4;^(k,a;) = Ak, 




(B5) 


sS(k,c.) = sW(k,a;) = 5k, 




(B6) 


4^i(k,c.) = s(^j(-k,-c.)=c,k+(|) j:^iiiiiiiu,-. 








T/(4) 2 T/(6) 2 
k,p,q,[k+p-q] k,p,q,[k+p-q] 




, (B7) 




-d; — En — Ea — -Ek+n-n + ^(^ CJ + -En + -^a + -Ek+n-a 


-2^] 



sS(k,c.) = sg(-k,-^) = C2k+(|)'E2^k^p^SCp-qSgn(7G) 



pq 



X V 



(4) 



,K 



(6) 



2{Ep + Eq + -Ek+p-q) 



[B8) 



k,p,q,[k+p-q] k,p,q,[k+p-q]^2 _ (^^ + Eq + E;k+p-q)2 ' 

where [k + p — q] is mapped to (k + p — q) in the first BZ by the reciprocal vector G. 
Feynman diagrams for the second-order self energies are shown in Fig. [TBr b)-(c). Above the 
coefficients Cik and C2k are 



Cik = ^^'(t?) J2^lll\-Glii2-l-k)xkXiX2 + ^i{2)xlx2 + ^i{2)xlxjx2 



12 

„2 I ,, /i\^2. 



+ 27i(k)xkXi + 7i(l)XkXi + 7i(2)x2 
2^Kn) 



^^ ^"^ )/^(l + 4)Ck, 



4V1 + C 



(B9) 



C2k = T^^lilTr) X]^i^2 37i(2- 1 -k)xiX2 + 37i(2- 1 -k)XkXiX2 -27i(l)xkXix: 



12 



27i(2)a;ka;2 - 7i(k)a;2 - 7i(k)xkX2 



1 / 2r] 



2V1 + C 



/k/rikCk, 



(BIO) 



with 



Ck 



2\2 



N 



J2 I'M [272(k) + 72(1) + 72(2) - 472(k +1-2) 



12 



xl 



+ 



72(2) -72(1 + 2 -k) {l + xixi)\. 



(Bll) 



33 



The divergent terms in Cik and C2k for k — )■ are canceled out by the second terms in 
Eqs.|B7]and|Bii^ 

Appendix C: Second order Coefficients and Vertex factors for the CAF phase 



The coefficients that appear in the second-order corrections in the Hamiltonian for the 
CAF-phase are: 

" ' E ^i'^2' I - 67i(2 - 1 - k)4a:;4 + 7;(2)xf 4 + I'A'^Xx'^ 

1/ C 



a 



Ik 



1 / 2 \ 2 

2HiV. 



12 



„'2^/2 / 
•^2 



+ 27i(k)x;,x'/ + 7i(l)4'^i + 7l(2)x^ 



4Vl + 2r7 



m+<)c'^ 



(CI) 



C^2k = \lk (^)' E ^'^2' [37U2 - 1 - k)x;4 + 37; (2 - 1 - k)x;,2x;4 - 2i,{l)x'^x\x' 
12 



2" ViV, 

27i(2)xkX2 - 7i(k)a;2^ - 7i(k)a;kxf 



1 



c 



2\l + 2-q 



^k^k^k) 



where 



and 



/,' 



1 + el 



2e 



1/2 



m. 



-sgn(7k^ 



k -• 



2e(, 



1/2 



^k'k) 



(C2) 



(C3) 



k -' 



a 



12 



27^(k) + 7^(1) + 7^(2) - 47^(k +1-2) 



+ 



(C4) 



7^(2)-7^(l + 2-k)](l + a;'M2)}. 

For the magnetization and spin-wave dispersion calculations only the vertex factors V'^^"* 
and V'^^^ are required which are: 



V, 



'(4) 



1234 



-7^(2 - A)x', - 7^(1 - A)x[x',x', - 7^(2 - 3)a;'3 - 7^1 - S)x[x',x', 



1 

+ 2 



7^(2) + l[{l)x[x', + 7; (3)44 + 4(4)44 
+ 7^(2 - 3 - 4)4^; + 7^(1 - 3 - 4)x[x'^x'3x'^ + 7^(3 - 2 - l)x[x'^ + 7^(4 - 2 - l)x[x'^ 
C 



1 + 277 



/1234 -^2 "I" sgn(7QJX]^X3X 



4 



(C5) 



K 



iS = 4(2 - 4)x'24 + 4(2 - 3)44 + 4(1 - 3)44 + 4(i - 4)44 



- - [-f[{2)xW3x'^ + 4(3)4 + 4(2 - 3 - 4)4 + 4(3 - 2 - l)x[xW4 

+ 4(1)444 + 4(4)4 + 4(1 - 3 - 4)4 + 4(4 - 2 - 1)444 
C 



l + 2r? 



4234 2:30;; + sgn(7^)x;a;2 



(C6) 
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with 

/i234 = ^ 72(l-3)+72(l-4)+72(2-3)+72(2-4)-7^(l)-7^(2)-7^(3)-7^(4) . (C7) 
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